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Abstract 


A technique for predicting the supersonic pitch damping of blunt axisymmetric bodies from static CFD 
data is presented. The contributions to static pitching moment due to forebody and aftbody pressure 
distributions are broken out and considered separately. The one-dimension pitching moment equation is 
cast to model the separate contributions from forebody and aftbody pressures with no traditional damping 
term included. The aftbody contribution to pitching moment is lagged by a phase angle of the natural 
oscillation period. This lag has been predicted by others using computational methods; it represents the 
time for aftbody wake structures to equilibrate while the body is oscillating. For a case with constant 
forebody and aftbody pitch-moment slopes, the characteristic equation of this formulation indicates that the 
lagged backshell moment adds a damping moment equivalent in form to a constant pitch damping term. 
The nonlinear moment equation is then linearized to yield an expression for pitch damping at a local angle- 
of-attack that is a function of the local contribution to pitching moment by the backshell and a lag time. 
CFD calculations of a capsule backshell’s contribution to the static pitching moment for a range of angles- 
of-attack is used to predict pitch damping coefficients. These predictions are compared with ballistic range 
data taken of the Mars Exploration Rover (MER) capsule and forced oscillation data of the Mars Viking 
capsule. The lag model appears to capture dynamic stability variation due to backshell geometry and Mach 
number. This lagged backshell pitching moment method can be used to make preliminary estimates of the 
dynamic stability of a blunt body traveling at supersonic speeds. 


Nomenclature 


a 

Speed of sound 

s 

Reference area 

B 

Linearized solution constant 

t 

Lag time 

Cm a 

Pitching moment coefficient 

t(f) 

Lag time 

Ctti^^bs) 

Backshell contribution to total C ma 

V 

Velocity 

Cm & ( F B) 

Forebody contribution to total C moi 

X 

Arbitrary equation variable 

Cm a 

Pitch damping coefficient (^gf) 

Greek 


Cm Q 

Linearized C m offset 

a 

Angle of attack 

d 

Diameter 

OL 0 

Linearization angle of attack 

E 

Linearized solution constant 

7 

Linearized harmonic oscillator offset 

f 

Frequency (cycles per second) 

e 

Angle of attack deviation from a Q 

Cy 

Moment of inertia about pitch axis 

A 

Characteristic equation roots 

Mbs 

Dimensional pitching moment (backshell) 


Harmonic oscillator damping coefficient 

Q 

Dynamic pressure 

p 

Flowfield density 
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Greek (cont’d) Subscripts 

4> Wake lag phase angle oo Freestream 

uj Harmonic oscillator natural frequency (rad/s) Overscript 

uj 0 Harmonic oscillator spring constant ' 1st time derivative 

Forebody spring constant “ 2nd time derivative 

CJ 2 Backshell spring constant 

Introduction 

The dynamic stability of blunt bodies is an important area of research for planetary entry mission 
design. As an entry capsule descends to low supersonic speeds (Mach_^_3.0 to 4.0 and below typically), base 
pressure contributions to the pitching moment become significant and can result in dynamic instabilities. 
These instabilities cause oscillations in angle-of-attack to increase in amplitude, potentially creating problems 
for parachute deployment and in extreme cases can cause the vehicle to tumble. It is therefore desirable to 
characterize the pitch damping of a blunt entry vehicle as well as possible in the low supersonic regime to 
mitigate any flight performance risks through entry trajectory modifications or stabilization devices like a 
supersonic drogue parachute. 

There are several test methods available to characterize the pitch damping of blunt vehicles. Forced 
oscillations tests were used extensively for the Mars Viking missions, 1 Apollo and other vehicles. Forced 
oscillation tests offer several attractive features. Highr rate, accurate damping measurements can be made, 
and tunnel conditions and freestream gases are both controllable. Sting effects, however, raise the questions 
of whether the data is valid for a free flying capsule as the wake structure is thought to be the dominant 
source of dynamic instability. Some studies 2,3 have shown that sting effects can be mitigated to some degree, 
but they remain a source of error. Ballistic range testing is another attractive option for blunt vehicles. This 
type of testing has the benefit of a free flying model with no sting effects. In some facilities, 4 the freestream 
pressure and gas (air, CO 2 , etc.) can be controlled as well. However, extracting dynamic derivatives via 
ballistic range testing is a complicated and less definitive process than reducing forced oscillation data. 
Pitch damping is extracted through parameter identification techniques. The aerodynamic coefficients are 
determined by fitting simulated trajectories through time, position, and attitude data collected at discrete 
points as the model flies down the range. The constants in assumed functional forms describing the aero 
coefficients are varied to minimize the error between the simulated trajectory and shadowgraph data points. 
In such tests, Mach number and oscillation amplitude are varying rapidly as the model decelerates. With 
these variables changing significantly during each shot, and considering the relatively sparse data (up to 50 
photo stations), pitch damping cannot be resolved to the precision that forced oscillation data can achieve. 
Both of these techniques are time consuming and expensive. Numerical determination of pitch damping 
by way of time-accurate computation fluid dynamics (CFD) is a relatively new option and to date has not 
been validated. The ability to quickly and accurately determine the general pitch damping characteristics 
of a given blunt body would be a great benefit to entry aerodynamicists. It would allow more effective 
preliminary mission planning without expensive experimental programs, and would help reduce the scope of 
such experimental programs when higher fidelity is required. While not validated, some significant advances 
have been made with CFD that have shed some light on the physics of blunt body dynamic stability. 

Teramoto et al 5 have shown evidence that the physics driving the dynamic instability of a blunt entry 
vehicle is a lag of the base pressure distribution behind its instantaneous attitude as the body undergoes 
oscillatory motion. They show computational results where the annular vortex in the wake flow behind an 
entry capsule takes a finite amount of time to reach a steady state structure when the body attitude is 
changed. For oscillatory motion, the wake structure lags by this amount of time behind the current body 
attitude. For the case they presented, this lag results in a dynamic instability (neglecting the forebody 
damping). This idea of a lagging base pressure distribution is taken from this computational work and 
applied to available static CFD solutions in an attempt to predict the dynamic stability characteristics of 
the Mars Exploration Rover and Viking entry capsules. It was hypothesized that static measurements or 
CFD predictions of backshell contributions to static stability might be an indicator of how dynamically 
stable or unstable a configuration might be. This work shows that lagging the backshell static moments does 
predict damping behavior that is consistent with experimental data. 
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Analysis 


Consider the blunt body in Figure 1 traveling at a constant velocity, V^. The body is free to oscillate 
and is statically stable. This configuration could represent a free-to-oscillate wind tunnel test where the 
tunnel velocity is constant and the model is constrained to rotate only about the sting axis. The pressure 
distribution on the forebody dominates the static aerodynamic forces and moments. The forebody also 
provides positive damping when oscillating, but the separated wake flow can introduce dynamic instabilities 
at small angles. As described earlier, it is believed that a lag of the backshell pressure distribution caused by 
the inertia of the separated wake structure is the mechanism that produces the instability. That phenomena 
is what will be modeled here. 



Figure 1. Blunt capsule traveling constant supersonic velocity, free to oscillate 


The equation of motion for this configuration is that for a simple harmonic oscillator and textbooks 
readily provide solutions for the angle of attack history of the body for constant coefficients. Other solutions 
exist for conditions where a body is free to heave due to lift and free to decelerate due to the drag force acting 
on the vehicle. 6 The special case of constant velocity with no heaving due to lift is used here to predict the 
pitch damping coefficient from static aerodynamics. A setup with more degrees of freedom would not change 
the contribution to amplitude growth from the pitch damping term, but additional degrees of freedom would 
add terms to the equations that are irrelevant to this analysis. The pitch damping predicted here should be 
relevant to other flight conditions with additional degrees of freedom. 


Linear Analysis 

The 1-D pitching moment equation for a blunt body constrained to constant velocity and a specified center 
of rotation (typically for a wind tunnel test, this is at or near an expected center of gravity) is given below. 
Summing the moments yields: 

!yy & ~ \pVlSdCm a A-d _ ^pV^SdC ma a = 0 (1) 

Which can be re-expressed in the form: 
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From the classic solution of a harmonic oscillator, substituting 
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a = e 


( 5 ) 


At 

into equation 2 yields the 2nd order characteristic equation 

(A 2 + 2£A + cj 2 ) e xt =0 => A 2 + 2£A + c j 2 = 0 (6) 

The roots of equation 6 provide the solution for equation 2 when substituted back into equation 5. Here, 
the forebody contribution to pitch damping is neglected by setting the usual damping term, £, to zero in 
Equation 2. This is not typically a valid assumption, but is done here to more clearly illustrate the effect 
of lagging the backshell contribution to the static pitching moment. A nonzero value for £ to account for 
the forebody damping will be incorporated into the solution developed in this section later in the paper to 
estimate the total damping of some example configurations. 

Now the pitching moment slope is separated into forebody and backshell contributions with the assump- 
tion that both are constant. 


Cm a — Cm (y{FB) + Cm (yiBS) (7) 

By imposing a lag in the backshell contribution to pitching moment as described by Teramoto et al, the 
moment equation can be cast in the following manner 


a(t) + l o\ot(t — t^) + l o\a(t) = 0 

Where 


( 8 ) 


and 


C, 




qSd 


( 9 ) 
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CO 1 


rn a(FB) 


qSd 


Again we have the exponential form of the solution given in Equation 5. 
introduce the lag time into the assumed form of the angle of attack 


_ tjL \ _ p \ (t-t#) _ e \t e -\t 


a(t — t^) = e 


Substituting Equations 5 and 11 into Equation 8 yields 


( 10 ) 

For the backshell term we 

(ii) 


(A 2 + u>le~ xt * + u>l) e xt = 0 =► A 2 + u 2 2 e~ xt * + w 2 = 0 (12) 

The exponent in the lagged pitching moment term in Equation 12 can be expressed as a truncated 
Maclaurin series. Recall the series expansion: 

e x = 1+X + — + — + ... (13) 

For small values of A the first two terms of the expansion are a reasonable approximation. Substituting 
into Equation 12 yields 


A 2 T (^ 2(1 — At^) T l j 2 — 0 

Equations 4, 7, 9 and 10 can be used to combine the UJ 2 and terms. Rearranging we have 

A 2 — LO^tcf) A + L d 2 = 0 (14) 

Equation 14 is now in the exact form of Equation 6. This shows that with the few approximations and 
assumptions made, a lag of the pressure distribution on the backshell of a blunt vehicle creates a damping 
behavior. Comparing the second terms in Equations 6 and 14, a relation for pitch damping can be determined. 
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L yy 


The pitch damping can then be determined from the base contribution to the static pitching moment 
and a lag time. 






OO t(j) 


(15) 


With some manipulation, we can get a more physical sense of the parameters influencing the damping 
effects due to the backshell. First, the lag time will be expressed as a phase angle or fraction of the period 
of one oscillation cycle. The roots of Equation 14 include the natural frequency of oscillation. These roots 
are: 


Ai ? 2 = — £ ± iuj 


where the natural frequency, cj, in rad/s is 


uo 


= V^o~C 2 = \ ~ 



C ma qSd 


L yy 
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L yy 


qSd 
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(16) 


The Cm a term dominates Equation 16 for our cases where is small and Cm a(BS) is small compared to the 
total pitching moment slope. Therefore, the oscillation frequency is roughly the same as the undamped value, 
uo ~ yj —CmaqSd/ lyy. Now, the lag time can be expressed as a fraction of a complete cycle of oscillation. 
Recall the conversion of frequency in radians per second to cycles per second: 


/ = 


UJ 
27 r 


t(f) — 
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(17) 


Substituting Equation 17 into Equation 15 yields a simple expression for the upper bound of the pitch 
damping due to the backshell pressure distribution lagging behind the capsule attitude: 
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V 7 - Cm a V PooSd 3 


(18) 


This expression requires only some reference parameters specific to the vehicle, static aerodynamic data 
and a phase angle to determine the pitch damping at a given atmospheric condition (freestream density). 
This equation for the pitch damping due to backshell effects, when substituted into Equation 3, yields the 
following expression for the moment due to damping (second term in Equation 2): 


MbS— damping 


= 2£b = -Q 


m a(BS) 


qSd 

C T ( 

^rria-Lyy 


(19) 


This equation shows that the damping due to the backshell is directly proportional to both the pitching 
moment slope of the backshell as well as the phase angle of the pressure distribution lag time. The damping 
is also related to the square root of the geometric parameters of the capsule (area, diameter) as well as the 
dynamic pressure at which the vehicle is flying and inversely proportional to the square root of the static 
stability of the forebody and moment of inertia. 

To predict the pitch damping for a particular configuration, static CFD solutions of a capsule along with 
its external dimensions and mass properties are all that are required except for one remaining item, the lag 
phase angle. Estimation of the lag time will be addressed below in the context of an example case. 
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Nonlinear C m (BS ) Analysis 

The analysis shown thus far assumes that C ma(BS) is constant. For the wake contribution to C m for a blunt 
body, this may only be appropriate at small angles-of-attack. Linearizing Equation 8 about an angle-of-attack 
will allow predictions of pitch damping to be made in nonlinear regions. 

First make the substitution: 


It follows that 


a = a Q + e 


( 20 ) 


a = 


e 


Where a Q is the angle-of-attack about which the moment equation is linearized and e is a small change in 
angle-of-attack from that point. For small deviations from the local angle-of-attack, the backshell pitching 
moment can be approximated by a linear curve: 


Cm{bs) ~ + C rYlo ~ + Cra 0 

The forebody contribution to pitching moment is still assumed to be linear. Substituting the linearized 
version of angle of attack and the linear approximation for the nonlinear backshell pitching moment into 
Equation 1, yields the following moment equation: 
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( 21 ) 


We again neglect the C m& term here and introduce the lag time to the backshell contribution to static 
pitching moment. Now, the lagged pitching moment equation (Equation 8) expressed in terms of e becomes 


C m lh , 
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which can also be expressed as 


where 


e(t) + oj\e{t) + l o\e(t - t^) + u 2 0 a 0 + 7 = 0 


7 = 


qSd 


Crr 


L yy 


( 22 ) 


Like the linear version, uq is the forebody pitching moment slope and remains constant in this analysis 
for all angles-of-attack. The form of Equation 22 looks similar in form to Equation 8 except for two constant 
terms which come from the linearization approximation. With a slightly different substitution for e we can 
recover the same characteristic equation. 

Let e be of the form 


e = e xt + B + E 

and the lagged backshell version of e is 


e(t - 4) = e xt e X C + B + 

Substituting e in Equations 23 and 24 into Equation 22 yields 

A 2 - | a„ t<p A + | a 0 +W„ \a Q B + | a „ 

Now the constants B and E can be set to the following values 


E 

E = w o \a„ a o + 7 
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B = —a, 


771 T G m 0 

-» 2 o~ C ma 


And Equation 25 reduces to 


A 2 — Cd2 (<^o) 2 ^0^ + ^o(^o) 2 — 0 (26) 

which is exactly the same as Equation 14 with the caveat that the constants are valid only near the angle 
of attack about which the equation was linearized. Again, the second term in Equation 25 can be used to 
determine an expression for the pitch damping due to a lagged base pressure distribution. 


n („ \ _ ^2(®'o) 2 VooIyyt'(j) _ ^rn a(bs) \ol q 

mM 0 ) ~ ” d 


(27) 


Substituting the simplified expression for the natural frequency into Equation 27 as was done with solution 
to the linear moment equation yields similar expression for the pitch damping due to the lagged backshell 
pitching moment slope, independent of freestream velocity. Just like the linear expression, the local nonlinear 
dynamic stability is a function of freestream density, vehicle parameters, the backshell and forebody pitching 
moment slopes and a phase angle: 


Crri^ (pt-o) 


2 a 


rn o i (BS) Wo 


2 L 


yy 




PooSd 3 


(28) 


Forebody damping 

For the complete description of the pitch damping of a blunt capsule, the forebody damping must be ac- 
counted for. This damping comes from differential changes to the pressure distribution over the forebody 
as the capsule oscillates (advancing portions see increased pressure and retreating portions see decreased 
pressure) and does not rely on any lag effects. Tobak and Wehrend 7 use a Newtonian method to estimate 
this damping for arbitrary shapes in the hypersonic flow regime. 

Adding the pitch damping term (recall £ as defined in Equation 3) back into Equation 25, we see that the 
forebody pitch damping can simply be added to the backshell contribution to arrive at the complete pitch 
damping of a capsule: 


A 2 + (2£ — CJ 2 \cto ^ 0 )^ + \a. 0 — 0 


Therefore 


n — C 

cx(total ) od. Newtonian) 


^CrriotiBS) Uo ^00 t(f) 


(29) 


Test Cases 

The mathmatical model presented will now be applied to a real case where static CFD calculations 
are used to “predict” the pitch damping behavior of two different blunt body shapes, the Viking and Mars 
Exploration Rover (MER) entry capsules. The results will be compared to pitch damping data obtained 
experimentally for these two shapes. Two different methods were used to obtain the pitch damping data, 
ballistic range for MER, and forced oscillation in a wind tunnel for Viking. Pitch damping data of both 
vehicles using both experimenal techniques would be more ideal for these comparisons, but the objectives of 
these comparisons are for general agreement only, not definitive validation of the approach. 

Figure 2 shows the dimensions of the Viking and MER outer mold lines (OML), each normalized by their 
maximum diameters. This figure shows the two backhells to be the primary difference between the outer 
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Viking 


MER 


Figure 2. Dimensions of Viking and MER capsules (nondimensionalized) 


mold lines. This difference and the existence of experimental data at similar conditions are the reason for 
their selection. 

/ 

MER Ballistic Range Data 

An example case will be used to perform a quantitative analysis. The example that will be used is a ballistic 
range shot of the MER entry capsule conducted at the Aeroballistic Research Facility at Eglin Air Force 
Base. 8,9 This case is used because it is relevant to this analysis, the dimensions and physical characteristics 
of the test models are well known, and pitch damping data were extracted from these tests with which 
comparisons may be made. The models tested were made of combinations of steel, tungsten and aluminum. 
The materials and proportions of each were varied to provide different center-of-gravity positions. They were 
fired down a 207 meter hallway instrumented with 50 shadowgraph stations that recorded their position and 
orientation. The time that each shadowgraph was taken was automatically recorded. This time and position 
information was then used in nonlinear parameter identification analysis to extract the static and dynamic 
aerodynamic coefficients. 

Details of the test facilities and the various shot conditions are discussed in in an earliar paper by the 
author. 9 The data reduction methodologies are described in more detail in a paper by Chapman et al. 10 The 
relevant dimensions, mass properties and flow conditions from those tests will be provided here along with 
results relating to the static CFD solutions calculated for this analysis. Table 1 lists the model dimensions, 
moment of inertia, freestream density, and approximate freestream speed of sound for the Eglin tests. The 
Eglin ballistic range facility is open to the atmosphere and therefore the speed of sound and air density varies 
slightly from day to day. These representative values are sufficient for this analysis. 


Table 1. Reference parameters for Eglin ballistic range and CFD predictions 


Parameter 

Ballistic Range Value 

d (m) 

0.07 

S ( m 2 ) 

0.00385 

>-'» H 

7? 

l 

to 

1.55E - 04 

Poo (kg/m 3 ) 

1.23 

a (m/s) 

348 
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Viking Forced Oscillation Data 

During the Mars Viking project in the early to mid 1970’s, extensive forced oscillation experimental data 
were obtained for the Viking aeroshell. 1 A portion of these data will be used to compare to results of the 
present analysis using CFD solutions of the Viking capsule. For comparisons to the CFD solutions, only the 
pitch damping results will be used. Parameters from the MER ballistic range testing will be used to calculate 
Cma values for comparison with these Viking data as well as the MER CFD predictions and ballistic range 
data. The use of one set of physical parameters and atmospheric conditions will simplify the comparisons 
of all these disparate data sources. The Viking project concluded that differences in the atmosphere in 
which the capsule is flying has little impact on the pitch damping of the vehicle. Viking also concluded 
that no corrections were required to account for damping data obtained from scaled vehicles to predict flight 
behavior. 11 These conclusions were subsequently validated by flight data from two successful entries and 
landings of Viking capsules as well as the by the use of Viking data to predict the Mars Pathfinder entry 
behavior. 12 Ideally, to compare the results of this pitch damping prediction method to experimental data 
of two different geometries, the validation data would be of equivalent fidelity. MER and Viking data were 
selected even though the pitch damping data derives from two different experimental techniques because the 
forebodies of those vehicles are nearly identical. Therefore any differences in damping should be due to the 
backshell contributions (and in this case, test technique). Therefore, for the level of accuracy sought for 
the present analysis, it will be assumed that pitch damping values determined by ballistic range or forced 
oscillation are of similar fidelity and insensitive to capsule size and atmospheric differences (air versus CO 2 ). 
Again, the goal of the present analysis is to determine order-of-magnitude pitch damping characteristics and 
general trends driven by backshell geometry and vehicle velocity. 

LAURA Calculations 

Three angle-of-attack sweeps were calculated using the Langley Aerothermodynamic Upwind Relaxation 
Algorithm (LAURA). 13 LAURA is a finite volume, shock capturing algorithm capable of calculating flows 
in chemical and thermal non-equilibrium. The solutions generated for this analysis are all laminar, perfect 
gas, thin-layer Navier stokes solutions. The molecular weight and ratio of specific heats used for the various 
solutions were appropriate to Mars entry (43.26 kg/kmol and approximately 1.3 respectively). The MER 
CFD cases were those used in the MER aerodynamic database 14 and the flow conditions taken from the 
reference Mars trajectory used for those cases. Additional angles-of-attack were run to better resolve the 
backshell pitchign moment curve. The Viking CFD cases were run at very similar conditions with a grid of the 
Viking OML with similar fidelity to the MER grids. These solutions were used as they were readily available 
and relevant to available experimental data. The intent of this first look at CFD data was to determine if 
the static pitching moments on representative vehicles were of the correct magnitudes and slopes such that 
a reasonable lag time could produce dynamic instabilities similar to those observed in experiement. 

Table 2 lists the solutions calculated for each a-sweep. The pitching moment coefficients due to the 
backshell only were extracted from the pressure distribution over the surface cells of the backshell from each 
of these cases. The pitching moment coefficient curves determined from these solutions are plotted in Figure 
3. Note that for each condition and shape there is an initial steep static stability near a = 0°. The pitching 
moment decreases with a rapidly and reaches a minimum before shallowing and returning to almost no 
contribution at large angles. The location and magnitude of the most negative C m ^ point is different for 
each configuration but the general character of the backshell moment contribution is similar. 


Table 2. Angle-of-attack Sweep CFD Solutions 


Case 

2° 

4° 

6° 

8° 

11° 

12° 14° 

16° 18° 21° 26° 

MER(M=2) 

X 

X 

X 

X 

X 

X 

X X X X 

MER(M=3) 

X 

X 

X 

X 

X 

X 

X 

VikingfM— 3) 

X 

X 

X 

X 


X X 

X 
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Figure 3. LAURA calculations of backshell contribution to pitching moment coefficient. 


Phase Angle Selection for Backshell Pitching Moment Lag Time 

The time lag of the base pressure distribution is a key component to this methodology. However, a rigorous 
basis upon which a lag time can be estimated does not exist. The proper understanding of any lag of the wake 
structure should be determined from experiment, looking for sensitivities to Reynolds number, oscillation 
frequency, Mach number and other parameters. Currently, there are some estimates available from unsteady 
CFD by Teramoto et al. ? The interpretation of the CFD data gives some theoretical understanding of 
what might drive the wake pressure lag. However, experimental work is necessary to validate and better 
understand the CFD predictions. 

The Teramoto calculations showed that a 15° phase angle corresponds to a lag time of 0.0021 seconds 
for a 100mm diameter blunt capsule oscillating at 20 Hz (sinusoidal forced oscillation) at Mach=1.3. To 
demonstrate the pitch damping predicted by wake lag in this paper, the Teramoto phase angle (</> = 15°) 
will be used. To justify its use, let us look at the oscillation frequencies imposed or observed in the Teramoto 
unsteady CFD predictions, ballistic range, and forced oscillation tests. The Strouhal number will be used 
for this comparison. This nondimensional parameter is often used when scaling models for dynamic tests. 
The Strouhal number, defined in Equation 30, can be thought of as the number of oscillations that occur for 
a given distance traveled through the atmosphere (or amount of freestream air passing over a wind tunnel 
model) . 


St = 


cod 

2 ^ 


(30) 


The Strouhal numbers at the conditions where the lag prediction will be compared to MER experi- 
mental data are very similar to the value specified in the Teramoto CFD predictions: Streramoto = 0.015, 
StMER,M = 2 = 0.016 and StMER,M=s = 0.011. The Strouhal number at the Mach 3.0 Viking forced oscilla- 
tion condition is roughly twice the Teramoto value: Sty iking, m = 3 = 0-032. It is not clear what parameters 
the phase lag is sensitive to or to what degree, so the Teramoto value of 15° will be used for all comparisons. 
It is reasonable to think that, all things being equal, a capsule oscillating more quickly (higher Strouhal 
number) would see the wake flow lag by a greater phase angle, but how phase angle varies with reduced 
frequency (and other parameters like Mach number) is beyond the scope of this paper. 
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With the selection of a phase angle for these comparisons, we return now to Equations 27 and 28. 
The local pitch damping coefficient can be calculated by multiplying the local pitching moment due to the 
backshell by single constant. For the MER ballistic rang conditions the expression for local pitch damping 
due to the lagged backshell static stability curve becomes: 


Cm«= — 


^m^BS) Uo Vootcf) 

d 


\<*0 / 2 Iyy 




PooSd 3 


• = -24.i • c. 


matBS) \Oto 


(31) 


Forebody Pitch Damping Contribution 

For the complete dynamic stability of the MER and Viking capsules the contribution from forebody damping 
is required. The modified Newtonian method is extended to supersonic conditions for an approximate 
estimate of the forebody contribution to pitch damping as no analytical or empirical techniques have been 
developed for blunt supersonic flows. For a C p ^ max of 1.87 and a center-of-gravity position of x cg /D = .2498, 
the pitch damping of the MER forebody is approximately: 


a 


^'ct ( M od. N ewto 


-.316 


Range of Applicability 

Before using the relations developed above to to generate pitch damping curves, it is important to estimate 
the valid range of angle-of-attack for which these relations apply. Figure 4 shows two Mach=3.0 LAURA 
CFD solutions used to generate the pitching moment curves of the MER capsule. The first plot, the MER 
capsule at 4° angle-of-attack, shows a fully separated wake region with the separation streamline extending 
far from the backshell. The second plot in Figure 4 shows the MER capsule at 18° angle-of-attack. This 
case shows the shear layer is very close to the backshell, on the windward side, if not completely attached. 
For such a case, the a large portion of the backshell is no longer within a fully separated wake. For this 
analysis, the claim is made that as the backshell pressure distribution becomes significantly influenced by 
attached flow and the wake ring-vortex structure breaks down at higher angles-of-attack, any lag in the base 
pressure has diminishing effect on the pitch damping of the vehicle. With the breakdown of the driving 
physical mechanism on the backshell, the forebody damping should dominate the vehicle damping behavior. 
The surface pressures on the backshell where the flow is attached will not lag like regions in separated wake 
flow. While the lag predictions of pitch damping will be made for all angle-of-attack solutions generated, 
this caveat should be noted when interpreting high angle-of-attack calculations. 

Results and Discussion 


Backshell Geometry Effects 

Cmoc predictions calculated using Equation 29 with the constant determined in Equation 31 are shown 
in Figures 5 and 6. The local slope at each data point was determined using an overlapping parabola 
interpolation scheme to provide smooth and continuous interpolated values between data points. A finite 
difference derivative was then calculated at each data point. Figure 5 shows calculated damping values for 
both the Viking and MER capsules at Mach 3. These cases are plotted against MER ballistic range and 
Viking forced oscillation data. Note that the C ma calculations were made at the center of gravity positions 
corresponding to the comparable data sources ( x cg /d = 0.27 and x cg /d = 0.214 for the MER and Viking 
cases). Overall, the agreement of the predictions with the experimental data is good. Both predictions 
indicate a region of decreased dynamic stability at lower angles-of-attack due to the lag of the backshell 
pressure distribution. In agreement with the experimental data, the lag calculations predict that the MER 
capsule is less dynamically stable than Viking at low angles. Quantitative conclusions about the higher 
angle-of-attack data are more difficult, but both predictions indicate increased dynamic stability at higher 
angles. Note that the ballistic range data is obtained by assuming a functional form and solving for the 
coefficients that model the behavior observed in the shadowgraphs. It is not surprising that the lag model 
shows local variation in the predicted damping curve that does not agree with the assumed form used in the 
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Figure 4. Mach = 3.0 LAURA solutions of the MER entry capsule 


ballistic range data reduction. In fact, one possible use of this new method would be to help with initial 
guesses of the functional form to used in ballistic range analysis. The Viking lag predictions show trends 
that generally follow the same form as the forced oscillation curve although there are obvious discrepancies. 
The lag calculations predict that the peak dynamic instability is not at 0 degrees. This agrees with the 
forced oscillation data and runs counter to the assumptions used in the ballistic range data reduction efforts. 
Considering the crudeness with which the lag time was determined and the potential breakdown of the lag 
analysis at higher angles-of-attack, the predictions are remarkably similar to corresponding experimental 
data. 

Mach Number Sensitivity 

Figure 6 shows two MER cases at Mach = 2 and 3 respectively, compared with ballistic range pitch damping 
curves. The lag predictions agree with the experimental data in one important area. At small angles, the lag 
calculations predict that the MER vehicle becomes more dynamically unstable with decreasing Mach number 
at lower angles-of-attack. The Mach 2 lag predictions show a much more significant peak instability near 5° 
than the Mach 3 peak. However, in terms of functional form, the Mach 2 calculations differ from the ballistic 
range data considerably more than the Mach 3 data. The point of greatest dynamic instability is again 
away from 0° and at larger angles there is an angle-of-attack range where the Mach 2 predictions indicate 
greater dynamic stability than the Mach 3 predictions. Again the method of selecting the lag time and the 
breakdown of the driving physics at high angles are likely sources of error. However, this prediction method’s 
ability to capture low angle-of-attack Mach sensitivity to pitch damping is an important finding. Comparing 
the static stability of blunt body backshells may be used as an effective screening method. Estimates of lag 
time might be imprecise, but when comparing backshell concepts for a candidate vehicle, identifying relative 
differences is quite useful even if experiments are still required to quantify the damping characteristics of 
the shape ultimately selected. The data evaluated here indicates that meaningful comparisons are possible. 
Further, the functional forms of these shapes may help with ballistic range testing, providing a more educated 
starting point for data reduction analysis. Conversely, a CFD-based functional form might be used in the 
data reduction, and the phase angle that best replicated the ballistic range data points would be extracted. 
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Figure 5. Cm^ predictions for MER (x cg /D = 0.27) and Viking ( x cg /D = 0.214) cases compared with ballistic 
range and forced oscillation data 



Figure 6. Cm & predictions for MER cases at Mach —2 and 3, compared with ballistic range data 
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Conclusions 


Linearized analysis of the constant velocity moment equation has provided a simple method of predicting 
the dynamic stability of blunt entry capsules. This method relies only on static measurement or CFD 
calculation of the contribution to the pitching moment slope by the vehicle backshell. Comparison with 
MER and Viking experimental data has shown that this analysis can predict dynamic stability trends due 
to backshell geometry as well as Mach number. Dynamic stability characteristics were obtained for two 
blunt body geometries that agreed well with damping characteristics measured experimentally. The only 
unknown parameter required to generate the lagged pit eh- damping predictions was the phase angle by which 
the backshell moments lagged. The phase angle was taken directly from computational work by Teramoto 
et al. The comparisons with experiment data suggests this approach can be used in screening candidate 
backshell shapes and may assist in ballistic range data reduction. 

A more rigorous method of determining the lag time, perhaps based on the momentum of the vortex 
structures in the wake flow, is a recommended next step. The static solutions indicate that the mechanisms 
assumed to drive the supersonic dynamic instabilities break down at higher angles of attack. Verifying that 
this breakdown occurs, and understanding that process could lead to improved results at higher angles-of- 
attack. An experimental investigation of the backshell pressure distributions and wake flowfields behind 
these types of blunt body flows should be done to validate the pressure lag and underlying mechanisms 
identified with CFD by Teramoto et al. Flow visualization, surface pressure and moment measurements 
would rill be helpful. Obtaining experimental data to better anchor this simple predictive model has the 
potential to more fully understand the physical mechanisms that cause blunt body dynamic instabilities, 
which to date have only been characterized with empirical data, supported with important, but limited, 
computational investigations. 
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